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It is by now well understood that quantum dissipative processes can be harnessed and turned into a 
resource for quantum-information processing tasks. In this paper we demonstrate yet another way in 
which this is true by providing a dissipation-assisted protocol for the simulation of general Markovian 
dynamics. More precisely, we show how a suitable coherent coupling of a quantum system to a set of 
Markovian dissipating qubits allows one to enact an effective Liouvillian generator of any Lindbladian 
form. This effective dynamical generator arises from high-order virtual-dissipative processes and 
governs the system dynamics exactly in the limit of infinitely fast dissipation. Applications to the 
simulation of collective decoherence are discussed as an illustration. 


Quantum decoherence and dissipation have been 
regarded until recently purely detrimental to the aim 
of quantum information processing (QIP) [I1I2]. In¬ 
teractions with the environment in fact inevitably 
lead to entanglement between the quantum comput¬ 
ing system and uncontrollable degrees of freedom. 
This unwanted entanglement in turn results in a 
system subdynamics that is in general incoherent 
and irreversible: unitarity is quickly lost and with it 
the quantum information processing advantages e.g., 
computational speed-ups, one was seeking for. This 
state of affairs triggered a spectacular theoretical ef¬ 
fort that led to the discovery of a host of techniques 
to tame decoherence [3] as quantum error correc¬ 
tion HIS], decoherence-free subspaces [MS] , noise¬ 
less subsystems [MIS] and holonomic quantum com¬ 
putation [l4l [T5] . 

It is therefore a conceptually remarkable shift the 
recent realization that by reservoir-engineering dis¬ 
sipation can be harnessed and turned into a use¬ 
ful practical resource for QIP (see m for an early 
pioneering insight). For example one can dissipa¬ 
tively achieve quantum state preparation [la US], 
quantum simulation m, holonomic quantum com¬ 
putation [20] and even universal computation |21j . 
Simulation of highly non-trivial properties of matter 
as topological order [52] and non-abelian synthetic 
gauge fields [23] can also be accomplished by dissi¬ 
pative means. Finally, all forms of QIP that encode 
information in the ground state of a time-dependent 
Hamiltonian, e.g., open system adiabatic quantum 
computation and quantum annealing, also benefit 
from dissipation and relaxation to negate thermally 
driven errors [241126] . 

In particular in |27j it has been shown that quan¬ 
tum information can be encoded in the set of steady 
states (SSS) of a sufficiently symmetric strongly 
dissipative system and manipulated coherently by 
an effective dissipation-projected Hamiltonian. The 
latter is of geometric nature and is robust against 


some types of Hamiltonian and dissipative pertur¬ 
bations [ 55 ]. The key idea of Ref. m is a simple 
one: once the system is prepared in the SSS the 
fast dissipative processes adiabatically decouple non 
steady-states away while at the same time strongly 
renormalize the system Hamiltonian in such a way 
that the SSS remains invariant under this projected 
dynamics. This phenomenon can be thought of as 
a sort of environment-induced quantum Zeno effect 
[51 [301 at the superoperator space level 1281 . 



FIG. 1. A quantum system S (blue ball) is coupled with 
coupling strengths Qi to M qubits (yellow balls). Each of 
these qubits is subject to amplitude damping with rates 
T~^. Proposition 2 shows that in the limit of small n the 
qubits can be adiabatically decoupled and the effective 
dynamics of S is described by M Lindblad operators of 
strength g^n 

In this paper we extend the ideas of [27] to higher 
order. In the case in which the dissipation-projected 
Hamiltonian is vanishing, higher order virtual dissi¬ 
pative processes give rise, in a suitable limit, to an 
effective Liouvillian generator that leaves the SSS in¬ 
variant. However, at variance with the case studied 
in [23 this effective generator is no longer Hamilto¬ 
nian: a slow irreversible process unfolds within the 
SSS. We will show how this mechanism can be ex¬ 
ploited to the end of the simulation of any Marko¬ 
vian dynamics. More precisely, we will show that by 
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suitably coupling a quantum system to a structured 
reservoir comprising multiple qubits undergoing fast 
amplitude damping one can implement an effective 
Liouvillian generator in general Lindblad form |31j . 
We will illustrate our results by analyzing the dissi¬ 
pative simulation of qubits subject to collective de¬ 
coherence. 

Preliminaries. - Let Ti., [dim('H) < oo] denote the 
Hilbert space of the system and L('H) the algebra 
of linear operators over it. A time-independent Li¬ 
ouvillian super-operator Cq acting on L(H) is given. 
The SSS of Cq consists of all the quantum states p 
[p G L(H), p > 0, Trp = 1) contained in the ker¬ 
nel Ker£o := {A1 / £o(-^) = 0} of Cq. We shall 
denote by T^o = ’^o (2o := 1 — ’^o) the spectral pro¬ 
jection over Ker£o (the complementary subspace of 
Ker£o-) As in m the Liouvillian £o is also assumed 
to be such that: a) (t > 0) defines a semi-group 
of trace-preserving positive maps with 116*"^° || < 1, 
b) The non-zero eigenvalues A^, [h> 0) of £o have 
negative real parts, i.e., the SSS is attractive. In this 
case Vq = limt_>oo 6*^° and Vq Cq = ^o'Pq = 0. We 
also denote by S := — lim^_,.o (2 — £o)~^Qo the re¬ 
duced resolvent of Lq at {z = 0) and by th := |j5||. 
The latter provides a natural time-scale associated 
with the relaxation processes described by Lq- The 
energy scale is of the order of the dissipative 
gap of Cq i.e., the smallest modulus of a non-zero 
eigenvalue of Cq. The dimensionless (and normal¬ 
ized) resolvent is defined by S := t^^S. We now add 
an Hamiltonian term K, := —i[K, •] where K = 
such that Lt = £o + A. We set K, = in 

such a way that A is dimensionless and ||A|| = 0(1). 
The time-scale T is our scaling parameter and has 
to be thought of as large or even infinite in the spirit 
of the adiabatic theorem. We first establish 
Proposition 1: If VoJCVo = 0, then for sufficiently 
large T one has that 

sup ||(e*^--eT^-)7^oll (1) 

tG[0, ST] V ^ 

where CeS '■= -VolCSiCPo, Ct '■= £o + ^ ^ 

0(1) > 0 and Cg = 0(1) depends on 9 and Cq. 

Proof. - Is provided in the Appendix. □ 

This result provides the starting point of this 
paper. In particular from Eq. Q it follows that 
limT-j-oo II— e®'^'=“')Po|| = 0- In words: if the 
system is prepared at time t = 0 inside the SSS and 
then evolves for finite fraction 6 of T, in the large T 
limit the time-evolution leaves the SSS invariant and 
it is governed by the effective (dimensionless) gen¬ 
erator £eff- [It is also sometimes convenient to in¬ 
troduce the effective dimensionful generator £eff := 


—PqJCSJCPo whose norm is |l£eff|| = 0{\\K\\^Tn). In 
terms of £eff the second term in the norm of Eq. Q 
reads e*’'''®.] 

Remarks: 0) The stronger the dissipation outside 
the SSS i.e., the shorter t^j, the weaker the effective 
one inside i) Since, by construction |l£eff|| = 0(1)) 
the (dimensionless) action associated to the effective 
propagator is 0{9) for T —oo. ii) The RHS of 
Eq. Q represents an error bound, if we fix it at e 
1, we see that one needs that T > CgT^. iii) If 
K, I— where ||Ai|| = 0(1) and VqICiVo 
0 then Q holds with £eff >-->• £eff -I- Vq^iVo and a 
different constant Cg =0(1) |32] 

The effective Liouvillian generator £eff is clearly 
reminiscent of the second-order effective Hamilto¬ 
nians routinely used e.g., in quantum optics, and 
obtained by some sort of adiabatic decoupling tech¬ 
nique |34j . However, this dynamics, at variance with 
that case as well as with the situation considered in 
m is not unitary but of general Liouvillian type. 
The key point is that this effective non-unitary dy¬ 
namics depends on K. and on its non-trivial inter¬ 
play with the bare dissipation generated by Cq . This 
opens up the possibility of using it to engineer dissi¬ 
pative systems with a desired Liouvillian generator. 

Universal Lindbladian simulation.- Let us con¬ 
sider a system S coupled to a system B via the gen¬ 
eral Hamiltonian 

M 

K = (2) 

i=l 

where the tensor ordering follows that of the total 
Hilbert space, H = Rs ® Rb and, without loss of 
generality, we assume b\ = Bi, {i = 1,..., M). We 
also assume that the dissipative term is of the form 
Cq = I 5 0 such that Cb(,Pq) = 0 where pq is by 
assumption the unique steady state of £ 3 . The SSS 
of Cq is given by all the states of the form p®pQ and 
it is isomorphic to the full-state space of S. In this 
case one has Vq{X) = TrB(A)(g)po and Vq1CVq{») = 
•] with ATeff = Tr^ (ATpo) 0 1 b[^- Let Sb 
be the the projected resolvent of Cb at 2 ; = 0 . 
Proposition 2: If ATeff = 0 then CeS = £gg^( 8 )lB, 

ip) ~ ~*[ILeff)P]+ ^Ty(LipAj —-{LjAi,p}) 
ij = l 

(3) 

where r := (r(^)+r(At)/2, = ^E^=i(r(^)- 

and = -Tr {Sb{B,pq)B,) . 

Proof. - Is provided in the Appendix. □ 

Notice that hI^ = Heff and that Eq. (j^ describes 
a truly Lindbladian dynamics iff T > 0. Our main 
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result now follows as a particular case of Proposition 
2 above. Let us consider a d-dimensional system S 
coupled to a system B comprising M qubits, by the 
Hamiltonian K = J2iLi9i ® + h.c.), where 

the Li’s are given operators acting on the system 
state-space only. Let us also suppose Cb = 
where each of the M qubits independently dissipates 
according to the local Liouvillian 

^i{p) = P^t - > P})- ( 4 ) 

The unique steady state oi Cb is po = |0)(0|® and 
since CT:{af po) = 0 (Vi) one has LTeff = 0. 

(S) 

Proposition 3: CeS = ^17 ® 1b where 

M 

4ff (p) = 4^5 ^,(L,pLI - -{lIl,, p}). (5) 

Proof. - To obtain Eq. ([^ from Eq. ([^, re-write the 
Bi,Lt in 1^ such that K = 9i (L|(8)cr“ -fh.c.). 

Remembering that po = |0)(0|®^, and Cb as in 
Eq. 0- we recover Eq. (|^ as required by direct 
evaluation of the matrix T™) in Prop. 2 □ 

Notice now that, in view of remark iii) after 
Prop. 2, one can add any Hamiltonian /Ci = T~^ICi 
[IIACill = 0(1)] acting on the system S only (^ 
VqICVo = V^fCiPo). This will result in £eff 
£eff -I- ICi- Therefore we see that Prop. 2 shows that 
in principle any Liouvillian in the Lindblad form 
|31| i.e., the most general generator of semi-groups 
of Markovian CP maps, can be obtained given the 
availability of M auxiliary qubits (one for each Lind¬ 
blad operator) subject to an amplitude damping 
channel and the ability to enact the Hamiltonian 
K. Dissipation turns into a resource that allows one 
to simulate a general Lindbladian evolution. 

We would like to make a couple of remarks: 1) One 
might think of obtaining the Lindbladian dynamics 
Eq. ([^ directly coupling the system S to some reser¬ 
voir with an interaction Hamiltonian of type ([^ and 
then using the standard Born Markov approxima¬ 
tion [34] . The point is that the latter involves un¬ 
controlled approximations (Markov) whereas Eq. (|T| ) 
has a uniform and controlled error 0(y'Tfl/r). 
This means that the effective dynamics of S be¬ 
comes exactly Lindbladian, with generator ([^, for 
T —oo. Of course this is true as long as the auxil¬ 
iary qubits are exactly described by the Lindbladian 
in Eq. 0 i.e., their genuine Markovianity is a key 
resource in our universal simulation protocol along 
with the ability of switching on the the Hamilto¬ 
nian in Eq. ([^. In view of physical applications, 
we stress that the effective dynamics in Eq. 0 still 


holds if the M qubits are replaced by M bosonic 
modes subject to amplitude damping i.e., the a~ in 
Eq. Q are replaced by annihilation operators Oi. 

In the next section we will discuss, for the sake of 
illustration of our general results, the simulation of 
different types of collective decoherence when S is 
itself a set of multiple qubits. 

Simulating Collective Amplitude Damping. - Here 
we use our general result Eq. ([^ to simulate a qubit 
subject to collective damping. This type of symmet¬ 
ric noise is interesting as it admits decoherence-free 
subspaces EHH] and can be used to dissipatively pre¬ 
pare entangled states. Let us consider a system of 
N qubits coupled to a bosonic mode e.g., N atoms 
coupled to a cavity EM mode, via a (collective) 
Jaynes-Cummings Hamiltonian K = g -\-S'^a) 
moreover we assume that the system dissipates ac¬ 
cording to the Liouvillian Cq = Is Cb where 
Cb{p) =—ioj [a^a, p\-\-Tf7{cLp7 — \{a)a, p}). Using 
Eq. 0 with Li = S~ one finds the effective genera¬ 
tor ^^s\p) = ‘ig'^TB (5'“pS'+ - 4{S'+S'“, p}) where 
p is just the qubit state as the generator is trivial in 
the bosonic degrees of freedom (frozen at |0).) 
Proposition 2 shows that one can consider for the 
auxiliary qubits a Liouvillian that is more general 
than Eq. Q (as long as its steady state is unique). 
We illustrate this fact by considering the thermaliza- 
tion of an auxiliary qubit at non-zero temperature. 
Namely, we add to Eq. Q an excitation Liouvillian, 
such that now (p) = {cr~ pa'^ — ^{a'^a~, p})-l- 

{a~^pa~ — ^{a~a'^, p}). By explicit computation 
of the T matrix in Prop. 2 one can check that the 
new effective generator (in the system S sector) is 

(p) = E P}) (6) 

Q —± 

where = 4g^ ^ numerical check of 

the validity of Eq. Q is shown in Fig. 

Note that one can add a Hamiltonian term of the 
form Hb = ajB<x^ in such a way that the unper¬ 
turbed Lindbladian becomes a thermalizing Davies- 
type generator, for which the temperature is fixed 
by t+/t_ = exp(/?BWB) [33]. In a similar fashion, 
using the remark iii) above, we can add a prop¬ 
erly rescaled Hamiltonian term ICi in such a way 
that £eff Lgfj = £eff -I- VqICiVo, and C'^q is again 
of Davies type. In other words, also the effective 
generator can be made thermal and so it defines a 
temperature according to ~ exp(/3effa;eft) for 

some energy scale Weff. From ^ it follows 

Tbs/Tb = oJesl^B, namely one has cooling or heat¬ 
ing according to whether ujb > Weff or vice-versa. 
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FIG. 2. (Color online) Distance from the exact evolution 
(Et ■= and effective one with Liouvillian (j^ , as 

a function of Ij^/T. N = 3, r+ = 2, t- = 1, and g — 
(where the relaxation time is tu = )■ 

The linear fit is obtained using the least squares fitting 
on all of the data points, and the norm is the maximum 
singular value of the maps realized as matrices. 


FIG. 3. (Color online) Distance from the exact evolution 
{St := e*'^^) and effective one with Liovillian Eq. Q, 
as a function of logjQ(f). = 1, r+ = 2, r_ = 1, and 
g = {trT)~^^‘^ {tr = )■ Note that for the dashed 

line we have extended t past T, purely for convenience. 
The norm is the maximum singular value of the maps 
realized as matrices. 


Simulating collective dephasing. - Following a sim¬ 
ilar set-up as the previous subsection but with a 
Hamiltonian of the form K = g ^ , the effec¬ 

tive generator becomes that of collective dephasing 
along the a;-direction 

(P) = - 1{S^S^,P})- (7) 

T+ -h T_ 2 

In Fig. 1^ we plot the distance between the ac¬ 
tual and the effective evolution as a function of t 
for different time-scales T. According to Eq. 0by 
changing T cT (c > 1), we expect the distance 
to fall by a factor of yA (cf. in Fig. the maximum 
error falls from the dash to solid line by a factor of 
~ a/TO)- In the limit of T —oo, the exact evolution 
becomes identical to that of the effective one for all 
times t (the actual evolution ‘adiabatically follows’ 
the effective one). 

Conclusions.- There is increasing evidence that 
dissipative and quantum incoherent processes can be 
used to enact quantum information processing prim¬ 
itives, see e.g. [TM25] . In this paper we have shown 
how a suitable coherent coupling between a quantum 
system S and an environment comprising multiple 
qubits subject to strong Markovian dissipation, can 
be used to simulate universal Lindbladian dynamics 
over S. More precisely, by using high-order virtual 
dissipative processes, one can build an effective Li¬ 
ouvillian generator in arbitrary Lindblad form |31j 
that governs the dynamics of S exactly in the limit 


of infinitely fast dissipation. This approach has to 
be contrasted with the standard one in which Lind¬ 
bladian evolution arises from a weak Hamiltonian 
coupling to a unitarily evolving environment and in¬ 
volves uncontrolled Markovian approximations e.g., 
Born-Markov factorization of the joint density ma¬ 
trix [33]. We illustrated our results by numerical 
simulations of concrete physical models. Our find¬ 
ings show that Markovianity itself can be seen as 
resource in that it allows for universal simulation of 
an important class of quantum irreversible processes. 

Acknowledgements. - This work was partially sup¬ 
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Note Added. After the completion of this work we 
became aware of Ref. |35] where a different approach 
to simulation of Markovian systems is pursued. 


[1] W. G. Unruh. Maintaining coherence in quantum 
computers. Phys. Rev. A, 51:992, 1995. 

[2] D. Aharonov and M. Ben-Or. Polynomial Simula¬ 
tions of Decohered Quantum Computers. In Pro¬ 
ceedings of FOGS, p. 46, Los Alamitos, CA, 1996. 
IEEE Comput. Soc. Press. 

[3] D.A. Lidar and T.A. Brun (eds.). Quantum Error 
Correction. Cambridge University Press, 2013. 

[4] P. W. Shor. Scheme for reducing decoherence 
in quantum computer memory. Phys. Rev. A, 
52:R2493, 1995. 

[5] D. Gottesman. Class of quantum error-correcting 













5 


codes saturating the quantum Hamming bound. 
Phys. Rev. A, 54:1862, 1996. 

[6] P. Zanardi and M. Rasetti. Noiseless quantum 
codes. Phys. Rev. Lett., 79:3306, 1997. 

[7] D. A. Lidar, I. L. Chuang, and K. B. Whaley. 
Decoherence-free subspaces for quantum computa¬ 
tion. Phys. Rev. Lett., 81:2594, 1998. 

[8] D. Kielpinski et al. Decoherence-Free Quantum 
Memory Using Trapped Ions. Science, 291:1013, 
2001 . 

[9] E. Knill, R. Lafiamme, and L. Viola. The¬ 
ory of quantum error correction for general noise. 
Phys. Rev. Lett., 84:2525, 2000. 

[10] P. Zanardi. Stabilizing quantum information. Phys. 
Rev. A, 63:012301, 2000. 

[11] J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley. 
Theory of decoherence-free fault-tolerant universal 
quantum computation. Phys. Rev. A, 63:042307, 
2001 . 

[12] P. Zanardi and S. Lloyd. Topological protection 
and quantum noiseless subsystems. Phys. Rev. Lett., 
90:067902, 2003. 

[13] L. Viola et al. Experimental realization of noise¬ 
less subsystems for quantum information process¬ 
ing. Science, 293:2059, 2001. 

[14] P. Zanardi and M. Rasetti. Holonomic quantum 
computation. Physics Lett. A, 264:94, 1999. 

[15] L. M. Duan, J. 1. Cirac, and P. Zoller. Geometric 
manipulation of trapped ions for quantum compu¬ 
tation. Science, 292:1695, 2001. 

[16] A. Beige, D. Braun, B. Tregenna, and P. L. Knight, 
Quantum computing using dissipation to remain in a 
decoherence-free subspace, Phys. Rev. Lett. 85, 1762 
( 2000 ). 

[17] B. Kraus et al. Preparation of entangled states 
by quantum Markov processes. Phys. Rev. A, 
78(4):042307-, 10 2008. 

[18] M. J. Kastoryano, F. Reiter, and A. S. Sprensen. 
Dissipative preparation of entanglement in optical 
cavities. Phys. Rev. Lett., 106:090502, 2011. 

[19] J.T. Barreiro et al. An open-system quantum sim¬ 
ulator with trapped ions. Nature, 470:486, 2011. 

[20] Victor V. Albert, Stefan Krastanov, Chao Shen, 
Ren-Bao Liu, Robert J. Schoelkopf, Mazyar 
Mirrahimi, Michel H. Devoret, Liang Jiang 
Holonomic quantum computing with cat-codes 
arXiv:1503.00194 

[21] F. Verstraete, M. M. Wolf, and J. 1. Cirac. Quantum 
computation and quantum-state engineering driven 
by dissipation. Nat. Phys., 5:633, 2009. 

[22] C.E. Bardyn et al. Topology by dissipation. New J. 
of Phys., 15:085001, 2013. 

[23] K. Stannigel et al. Constrained dynamics via the 
Zeno effect in quantum simulation: Implementing 
non-abelian lattice gauge theories with cold atoms. 
Phys. Rev. Lett, 112:120406, 2014. 

[24] A.M. Childs, E. Farhi, and J. Preskill. Robustness 
of adiabatic quantum computation. Phys. Rev. A, 
65:012322, 2001. 

[25] M. S. Sarandy and D. A. Lidar. Adiabatic quan¬ 
tum computation in open systems. Phys. Rev. Lett., 


95:250503, 2005. 

[26] K. L Pudenz, T. Albash, and D. A Lidar. Error- 
corrected quantum annealing with hundreds of 
qubits. Nat. Commun., 5:02, 2014. 

[27] P. Zanardi and L. Campos Venuti. Coherent quan¬ 
tum dynamics in steady-state manifolds of strongly 
dissipative systems. Phys. Rev. Lett. 113, 240406 
(2014). 

[28] P. Zanardi, L. Campos Venuti, Geometry, robust¬ 
ness, and emerging unitarity in dissipation-projected 
dynamics, Phys. Rev. A 91, 052324( 2015) 

[29] P. Facchi and S. Pascazio. Quantum Zeno sub¬ 
spaces. Phys. Rev. Lett., 89:080401, 2002. 

[30] G. A. Paz-Silva, A. T. Rezakhani, J. M. Dominy, 
and D. A. Lidar. Zeno effect for quantum compu¬ 
tation and control. Phys. Rev. Lett., 108:080501, 
2012 . 

[31] G. Lindblad, On the generators of quantum dy¬ 
namical semigroups, Commun. Math. Phys. 48, 119 
(1976) 

[32] From the Proof of Prop. 1 reported in the Appendix 
it easily follows that only the term T~^T’oK-iPo is 
0(1) whereas all the other terms involving KA are 
at most 0{1/VT). 

[33] R. Alicki, Quantum Dynamical Semigroups and Ap¬ 
plications, (2007). Springer Science & Business Me¬ 
dia. 

[34] Gardiner, P. Zoller Quantum Noise, Springer 

[35] Ryan Sweke, Ilya Sinayskiy, Denis Bernard, 
Francesco Petruccione, Universal simula¬ 
tion of Markovian open quantum systems, 
arXiv: 1503.05028 (2015). 



6 


Proof of Proposition 1 

Here we provide a Proof of Eq. 0 and an asymp¬ 
totic (large T) estimate of the constant Cg. Let V 
be the spectral projection of Ct = Cq + 1C, (||/C|| = 
0{^) associated with the zero eigenvalue of Cq. 
Since, because of the Lindblad structure, there is 
no nilpotent term associated with the zero eigen¬ 
value, the perturbation theory reads, as shown in In 
T. Kato, Perturbation theory for Linear operators, 
for small ||/C||, i.e., large T, 

p - Po = -Po/C5 - 5/CPo + 0{\\ICf), 

VCtV = VolCVo - VolCSICVo - VqICVoICS 

-SlCVolCVo + OiWICf) (8) 

From the first equation it now follows (for sufficiently 
large T) 


definition £eff := —VqICSICVo for the dimensionful 
effective generator from the last equation in ^ it 
follows 

||P£T^-£eff|| =0(||/cf), (10) 

whence (for small ||Af||) ||P£t 7^|| < Cali^Ceffll- Since 

— gtvCT'P'p write 

{e*^^ - e*^^^^)Po = -(e*^-^ - e*^^^^)(P - Pq) 

gt£T _ gtVCTV ^ ('gtCT _ 

+ ( 11 ) 

Using lie''-— e^+^ll < with X := t£eff 

and Y = tlVCr'P — £eff) and the bounds above it 
also follows that, for 0 < t < 9T, 

< C^t||/C||3 (12) 


\\V-ro\\=0{Tn\\IC\\)<C[TR\\IC\\, (9) 

where is a suitable constant (notice that ||Po|| = 
1). On the other hand, using Pq/CPo = 0 and the 


where is a constant of that, for dimensional rea¬ 
sons, is O(t^) i.e., C 2 < C' 2 T^. From (11) using 
||e‘’^^|| = 1, and standard operator norm inequali¬ 
ties one finds 


e* :=||(e‘^^ - e*^»«)Po|| < ||(e‘^=« - e*^^^^)Po|| + (||e*^^|| + ||e*^^^^||)||P - Po|| 
< + {1 + e‘ll^^'^^ll)||P - Po|| 


(13) 


Notice that et is the quantity showing up in the LHS 
of ([^ namely is the quantity whose upper bound 
over [0, T] we desire to show is 0{^/trJT). Now 
using the bounds ([^, (10),(12) and 0 < t < dT, {9 > 
0 ) one finds 


for t G [0, 9T] one obtains sup^ et < \/^ (C'i-I-C '2 9). 
Setting Cg \= C\P 0^9 completes the Proof of ([^. 


Proof of Proposition 2 


et<TR\\IC\\{Ci+C2tTR\\ICf), 


(14) 


where Ci > C[{\ + ). By moving to 

the dimensionless Hamiltonian such that ||/C|| = 
{TRT)-^l^iC\\ the inequality (|l^ becomes et < 


y/^{Ci + ^C 2 )- Notice that the requirement of | 
being sufficiently small used repeatedly in the above 
translate now in the “adiabatic criterion” of T being 
sufficiently large. Finally by taking the supremum 


We directly compute the second order effective 
generator '■= —VqICSICVo by acting on some 
state X, such that Vo{X) = p ® pq. One has 
SICVoiX) = -ij:f^^iL,p®SBiB,po) - 
pLi 0 SBipoBi)), where we have introduced nota¬ 
tion S = Is (E> Sb, which only acts non-trivially on 
system B for this set-up, i.e. Sb = — /q°° (fol¬ 
lowing from S = — e*^°Qo, see [17]). Acting 

with —VqIC on this we can see that: 


M 


M 

C,s{x) = { V {L,pL, - L,L,p) + V {L,pL, - pL,L,)} 0 po. 


hi=i 


hi=i 


(15) 


where P-^^ = —YT{SB{BiPo)Bj), and P^^^ = can rewrite the system S'part of these equations in 

—Tr {SB(pQBi)Bj). In passing we notice that one 
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a more familiar form, using that without loss of gen¬ 
erality Li = l\ and Bi = b]. Just observe that 
since iSs is a Hermitian-preserving map, we have 
r("4)* = Eq. ^ then follows as required. 


The dissipation-projection hierarchy 

Before concluding, we would like to show how the 
construction leading to the effective dynamics Q 
can in principle be iterated over a sequence of ex¬ 
ponentially longer time-scales. Let us set the er¬ 
ror parameter e := |jiL||T/j ^ 1. The effective re¬ 
laxation time of Q can be roughly estimated as 

4^^ IlCft II ^ » Ti^, where 

the last inequality stems from the condition e <§; 1. 

Suppose that the dynamics generated by CeS ad¬ 
mits itself a high-dimensional SSS (let denote 
the associated projection) and that one can switch 


on an extra Hamiltonian Ki such that = 

e ^ 1. One can now apply the projection the¬ 
orem 0 to £eff and Ki and argue that the ef¬ 
fective dynamics in the SSS of £eff is ruled by 
/Cgg := ICiVq^\ If even this effective Hamilto¬ 
nian vanishes then one can iterate the projection pro¬ 
cedure assuming as a starting state-space the SSS 
of £eff- In general, if at the n-th level one finds 
~ H 4"^ i® high-dimensional then 
one can move to the next level where = 

JCn'PQ^'^ and a new Hamiltonian iL„+i 
such that = e <C 1 is introduced. Rea¬ 

soning as in the above one can show that at each iter¬ 
ation the relaxation time scale (Hamiltonian norm) 
gets stretched (compressed) by a factor e~^ (e^). 
From the point of view of potential applications the 
interest in exploring this projection hierarchy rests 
on the possibility that the first non-vanishing ef¬ 
fective Hamiltonian has some desired property e.g., 
higher-locality |27) . 


